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ABSTRACT. After two decades of direct dynamical simulation of large— scale 
structure in the universe, it is safe to say the subject is now "mature". Still, 
there are parts of the problem that are less well developed than others. In 
general, the collisionless dynamics of the dark matter component is better 
understood than the collisional gas dynamics of the baryonic component. In 
situations where the gas dynamics is relatively simple, such as the Lyman— a 
forest and the intracluster medium in X— ray clusters, our ability to reproduce 
observational data has evolved rapidly, and the interpretive and predictive 
power of such experiments should now be taken seriously. A comparison of 
twelve gas dynamic codes to the problem of forming a single X-ray cluster 
shows that numerical inaccuracies are modest (typically below ten percent), 
leaving missing physics as the main source for large systematic differences 
between theory and observation. Galaxy formation, being more complex, is 
farther behind in its development, but simulations capable of resolving the 
morphological range of the Hubble sequence in cosmologically interesting vol- 
umes may be just around the corner. 



1 Introduction 

Building a universe in the labo- 
ratory is a daunting, if perversely 
inviting (think of the funding lev- 
els!), proposition. Fortunately for as- 
trophysicists, there are interesting 
physical questions about the universe 
that can be addressed without the 
creation of a bona fide mock-up — 
a reasonable facsimile will do. How 
"reasonable" the facsimile needs to 
be is dependent on the level of detail 
of the questions being posed and the 
complexity of the physics driving the 
processes at hand. 

In this contribution, I will briefly 
review the status of attempts at 
producing virtual facsimiles of our 
real universe through direct numer- 
ical simulation. In particular, I'll 
focus on the formation of large- 
scale structure (LSS), where "large" 
means galactic scales and upward. 
Sub-galactic scale evolution, notably 
the Lyman-a forest, is reviewed by 
Weinberg etal. in this volume. The 
reader will note a general theme 
which is both historical and epis- 
temological; "simpler" parts of the 
problem which have been worked on 
for the longest time are better under- 
stood than the more difficult aspects 
which are only now receiving careful 
attention. 



2 N-body Models of LSS 

On it's largest scales, structure 
in the matter component of the 
universe is driven by gravity. The 
first attempts at direct N-body mod- 
eling of the gravitational instabil- 
ity process in a volume of cosmic 
scale (Aarseth, Turner & Gott 1979) 
revealed morphological features — 
clusters, walls and voids — charac- 
teristic of the observed large-scale 
galaxy distribution. Future experi- 
ments enlarged the simulated vol- 
umes and increased the resolving 
power, leading to the current view 
of an evolving "cosmic web" of dark 
matter as seen in the images from 
the set of 256 3 particle Virgo sim- 
ulations (Jenkins etal. 1998) shown 
many times at this meeting. 

Since 1970, advances in comput- 
ing technology — faster processors, 
more memory, better algorithms — 
have allowed the typical number of 
particles in cosmological simulations 
to increase by a factor of 100 ev- 
ery decade. The increased dynamic 
range has been exploited in two or- 
thogonal directions which can be 
roughly characterised as increased 
volume at fixed resolution (minimum 
mass/length scales) or increased re- 
solving power (smaller minimum 
mass/length scales) at fixed simu- 
lated volume. 
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Table 1. Hubble Volume Models. Periodic cubes of side L are simulated using 
CDM fluctuation spectra with shape parameter T = 0.21 and normalizations 
erg listed. The redshifts correspond to look-back epochs along the box diagonal 
Z V3L' ^ ox s ^ e ^ en 9^ 1 Z L an d half-length z L /2 f or the light-cone datasets. 



Model 




Q A 




L 


V3£ 




ZL/2 


rCDM 
LCDM 


1.0 
0.3 


0.0 
0.7 


0.60 
0.90 


2000 
3000 


4.6 
4.4 


1.3 
1.5 


0.45 
0.58 



2 . 1 The Hubble Volume Project 

The natural progression toward 
simulating ever larger volumes has 
a practical endpoint — the entire 
visible universe. Though past ex- 
perimenters could have arbitrarily 
enlarged their volumes to encom- 
pass the Hubble length cH — 1 = 
3000 h" 1 Mpc, the lack of ability to 
resolve structure out of the linear 
regime of fluctuations left the exer- 
cise empty. Why simulate an analytic 
system? 

Parallel computers, with their 
large pools of multi— processor mem- 
ory capable of supporting very large 
N calculations, have changed the 
situation dramatically. Using 32- 
bit numerics, the phase space posi- 
tion plus an index (useful in post- 
processing) for 10 9 particles requires 
28 Gb of memory. A parallel machine 
with 512 nodes each with 128 Mb of 
memory offers 65 Gb, enough for 10 9 
particles with space left over for ad- 
ditional large arrays needed for book- 
keeping and gravitational potential 
calculations. 

What will a billionfl particles do 
for you? Suppose we ask that the 
particle mass be sufficiently small 
that the Coma cluster (mass ~ 
1.1 X 10 15 /i- 1 M Q ) be resolved by 
500 particles. Then the mass asso- 
ciated 10 9 particles, Mtot = 2.2 X 
10 21 /i- 1 M Q , would fill cubes of 
length 2000 h' 1 Mpc for Q m = 1 and 
3000 h' 1 Mpc for U m = 0.3. Since 
the diagonals of these cubes contain 
independent information on a length 
scale exceeding the Hubble length, 
they have a legitimate claim to be 
considered simulations of the "Hub- 
ble Volume" . 

Simulations using a billion par- 



ticles in such Gpc sized volumes 
have now been performed by the 
Virgo Consortiumn using a message- 
passing version of the Hydra N-body 
code (MacFarland et al. 1998) run on 
a 784 processor SGI/Cray T3E at the 
Rechenzentrum Garching. The runs 
each took roughly 5 days of cpu time 
on 512 processors (equivalent to 7 
years of single— cpu computing) and 
each generated ~ 200Gb of output 
data. 

With such a large simulated vol- 
ume, the traditional method of 
recording images of the mass dis- 
tribution at fixed intervals of world 
time is of little value for address- 
ing questions of observational signif- 
icance. Observational data comes to 
us from our past light— cone, so gen- 
erating output along the past light- 
cone of artificial "observers" in the 
simulated volume is a natural solu- 
tion. The Hydra code was modified 
to propagate output filters through 
the volume, creating "surveys" with 
different amounts of "sky coverage" 
to the depths listed in Table 1. 

Details of the geometries of the 
light cone output is available on the 
web[j], and a gif image of the "Hub- 
ble tie", a 40 /i — 1 Mpc thick slice 
along the diagonal wedge survey of 
the tCDM model can also be found 
there. This image, displayed in large 
format in the foyer of the MPA 
building during the meeting, reveals 
clearly the emergence of large— scale 
structure in the universe from z ~ 5 
to the present. 

Figure [l] shows the number of 
clusters found in each of the simu- 



1 pardon the Franco- American 
lingo. 



see http://star-www.dur.ac.uk/ 
Trazcrp/ virgo/ virgo. html 

3 http:/ /www. physics. lsa.umich.edu/ 
hubble-volume. Add /lightcones.htm 
to this address to get directly to the 
Hubble tie image. 
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Figure 1. Cumulative number of 
clusters identified in the Hubble Vol- 
ume models at the final epoch. 

latcd models at the present epoch. 
The shear number of collapsed ob- 
jects — roughly 3000 clusters at 
least as massive as Coma and nearly 
one million clusters resolved by 
more than 32 particles — allows 
detailed investigation of statistical 
questions such as the behavior of 
the cluster-cluster correlation func- 
tion with cluster richness (Colberg 
etal. 1998). 

Very rare events can also be stud- 
ied. In the tCDM model at the 
present epoch, the largest group 
identified with a standard percola- 
tion technique lies in a region where 
there are four clusters each more 
massive than Coma lying within a 
cube of side 20/i _1 Mpc! No obvi- 
ous counterpart in our local universe 
comes to mind (the Great Attractor 
region may come close), but future 
deep X-ray, Sunyaev-Zel'dovich and 
gravitational lensing studies should 
discover such rarities. 

A view of clusters resolved by 
at least 32 particles (M > 7.0 X 
10 13 /i -1 M Q ) and located in the 10 X 
10 square degree light— cone survey 
along the cube diagonal of the tCDM 
run is shown in Figure Q The plot 
shows positions in the background 
Robertson- Walker metric (the space 
of the computation) and the linear 
extent of the image is 2400 
Mpc. All the spatial information in 
the image is formally independent; 
no periodic replications are used to 
generate the map. 

A total of 1313 clusters exist in 
the field, the furthest located at red- 




1.5 



Figure 2. Positions of clusters in 
the 100 sq deg diagonal wedge light- 
cone survey from the tCDM model. 
Numbers indicate redshift. 

shift 2 = 1.68. The emergence of clus- 
ters at z < 1 is evident from the en- 
hancement in number density toward 
the vertex. Hints of the filamentary 
large-scale matter network are visi- 
ble at moderate redshift, but at high 
redshift the thickness of the wedge 
approaches 200 h~ 1 Mpc and projec- 
tion effects blend together the coher- 
ent features which have characteris- 
tic length scale ~ 40 h~ 1 Mpc. 

Basic observables of this cluster 
population are displayed in Figure [?|. 
The panels show the line-of-sight ve- 
locity of the cluster center of mass, 
the one-dimensional velocity disper- 
sion o~\r> and the cluster mass identi- 
fied using a friends— of— friends perco- 
lation algorithm with linking length 
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Figure 3. Characteristics of the 1313 
clusters displayed in Figure The 
minimum mass of 7 X 10 13 /i _1 M0 
corresponds to 32 particles. See the 
text for discussion. 



0.2 times the mean interparticle sep- 
aration in the comoving metric. 

A surprising result from this fig- 
ure is the fact that the cluster with 
the third highest velocity dispersion 
in the sample is also the most dis- 
tant object at z = 1.68. The one- 
dimensional velocity dispersion of 
a\D = 1072 km s — 1 translates into 
an intracluster medium temperature 
kT = 6.3 keV, assuming a ratio of 
specific energies o~^ D / (kT / fim p ) = 
1.17, (Frenk etal. 1998). The ther- 
mal Sunyaev— Zel'dovich (SZ) effect 
from such a cluster should be readily 
observable (see §3 below). The exis- 
tence of such a potential well at this 
redshift provides hope that such tar- 
gets will begin to appear in deep ob- 
servational searches for clusters using 
SZ, weak gravitational lensing, and 
even direct X-ray detection (which 
must fight (1 + z) 4 surface bright- 
ness dimming). A few examples of 
clusters with z > 1 have emerged, 
including a system discovered op- 
tically, CIG J0848+4453, lying at 
z = 1.273 with a velocity dispersion 
cr 1D = 700 ± 180 km s" 1 based on 8 



galaxy spectra (Stanford etal. 1997). 
Deep ROSAT HRI images of regions 
centered on high redshift radio galax- 
ies have also provided distant clus- 
ter candidates, among them 3C294 
at 2 = 1.786 (Dickinson, private com- 
munication). 

The upper panel in Figure jj| shows 
the line-of-sight (los) mean peculiar 
velocity v z of the clusters as a func- 
tion of redshift. This quantity is ob- 
servable via the kinematic SZ effect, 
which scales linearly with v z , but the 
expected amplitudes are small. Even 
for a cluster with v z = 1000 km s , 
the kinetic amplitude is smaller than 
the thermal effect by about a fact' 
10 (see Birkinshaw 1998). Figure 
displays a decrease in the character- 
istic los velocity with redshift, as ex- 
pected from linear theory, indicat- 
ing the kinematic SZ effect as a poor 
choice for detecting high-z clusters. 

Of course, the fact that these pro- 
cesses locally distort the cosmic mi- 
crowave background means that the 
CMB is the ultimate search engine 
for high redshift clusters. The future 
MAP and Planck missions, which 
will provide nearly all-sky coverage, 
will see bright clusters as foreground 
hot and cold spots, depending on 
wavelength, subtending roughly ar- 
cmin angular sizes. The lightcone 
data from the Hubble Volume will 
be used as a template for such ef- 
fects, with gas physics being intro- 
duced "by hand" into the potential 
wells in a model dependent manner 
(Kudlicki etal. , in preparation). 

The existence of a hot cluster 
at high redshift is assisted by the 
fact that the virial temperature at 
a fixed mass scales as T oc M/r oc 
p 1 / 3 (z) oc 1 + z (Kaiser 1986). 
This relation is shown in the mid- 
dle panel of Figure ^ for clusters 
at the 32 particle limit assuming 
a\ D = 1000(M/10 15 /i- 1 M Q ) 2 / 3 (l + 
z) (e.g., Bryan & Norman 1998). 
Clusters at the low mass cutoff fol- 
low this mean relation well (with an- 
ticipated ~ 10% scatter). 

The z = 1.68 cluster actually has 
mass a factor two above the 32 parti- 
cle limit. Is a 1.4 X 10 14 /i _1 M clus- 
ter expected in a survey of this size 
at this redshift? From Figure ^, it is 
evident that this object is a rare oc- 
currence, but it does not seem patho- 
logically so. The mass spectrum at all 
redshifts has a naturally ragged high 
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mass envelope, and the z = 1.68 ob- 
ject appears a reasonable extension 
of the trend exhibited at lower z. 

More quantitatively, one can use 
the traditional link between linear 
fluctuation amplitude and collapse 
epoch 6(z) = 1.68(1 + z) along with 
the employed rCDM normalization 
of the initial Gaussian fluctuations to 
estimate that this object corresponds 
to roughly a 5a upward density fluc- 
tuation. Now 5cr is certainly rare — 
the tails contain p = 2erfc(5) = 
6 x 10 — 7 of the Gaussian pdf — 
but countering this rarity is the large 
number of independent samples of 64 
particle masses available in the Hub- 
ble Volume simulation 7V = 10 9 /64 = 
1.5 X 10 7 . The product Np implies 
there should be about 10 such clus- 
ters in the entire simulation volume. 
Since the survey volume between z = 
1.5 and z = 2 in Figure |^ represents 
about 1% of the whole, it's a bit 
lucky (1 in ~ 10) but not crazy that 
one such cluster appeared in the sur- 
vey volume. 

2.2 Halo Structure 

One can use the increased dy- 
namic range from parallel comput- 
ers to probe the internal structure 
of individual objects such as clusters 
to higher density contrasts (Moore 
etal. 1998; Dubinski 1998). Com- 
pared to using the increased parti- 
cle number to "buy volume", this 
is a technically more demanding 
task. Since the gravitational dynam- 
ical time scales as (Gp)~ 1 / 2 , resolv- 
ing structure to higher densities re- 
quires integrating particle orbits us- 
ing shorter timesteps and over more 
dynamical times. 

Moore etal. (1998) break new 
ground with simulations in which a 
single cluster is resolved with over 
3 million particles interior to the 
virial radius (conventionally defined 
as r200i the radius of the sphere 
within which the mean interior mass 
density is 200 times the critical den- 
sity). With force resolution of 5 and 
10 kpc in two separate runs, they 
approach three orders of magnitude 
dynamic range in linear scale within 
the virial radius of 2 Mpc for their 
simulated "Virgo" replica (the name 
refers to the observed cluster, not 
the consortium). They find the in- 
ner portion of the density profile to 



approach a power law p(r) ex. r~ 1A , 
steeper than the asymptotic slope 
of —1 predicted by the model of 
Navarro, Frenk & White (1997 and 
references therein; hereafter NFW) 

4d = b. . (i) 

p c cx(l + ex) 2 

Here x = r/r2oo ls the scaled radius, 
c is the concentration parameter — 
the sole fitting parameter — and S c = 
(200/3)e 3 /[ln(l + c) - c/(l + c)] en- 
sures consistent normalization of the 
profile at x = 1. 

Moore etal. find that the NFW 
form does not provide an acceptable 
fit across the entire range of densi- 
ties probed by their experiment; for 
any value of the concentration pa- 
rameter, there remains structure in 
the residuals about the fit. The inte- 
rior portion of the profile appears to 
be driving the conflict. Fits to lower 
resolution runs are acceptable, but 
the inner part of the profile appears 
unwilling to cooperate and roll over 
to the r— 1 form expected from equa- 
tion^. 

Population aspects of The NFW 
profile were presented at this meet- 
ing by Jing (these proceedings) who 
examined a large body of halos from 
a set of 16.8 million particle runs 
performed on a parallel supercom- 
puter with 16 fast vector processors. 
Jing examines complete samples of 
halos within the simulations and gen- 
erates impressive statistics. A single 
run contains 300 to 500 halos re- 
solved by more than 10,000 particles. 
He finds that a significant fraction of 
objects have density profiles that are 
poorly fit by equation(Q). 

At first clance, this result is not 
surprising, since the studies by NFW 
selected against halos which were dy- 
namically active while Jing's study 
includes all objects. However, what 
is surprising is the size of the sam- 
ple which is not well fit by the NFW 
form; over a third of the clusters in 
the sample exhibit maximum frac- 
tional deviations of at least 35% from 
the NFW form (using 10 bins per 
decade in radius). The best fits for 
this fraction yield a distribution of 
concentration parameters quite dif- 
ferent from the remaining clusters in 
the sample — the peak concentration 
is lower by a factor two and the width 
is broader by a similar factor. 
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What does this all mean for the 
NFW profile? An obvious lesson to 
learn from these new studies is that 
one must exercise care in application 
of the density profile in equation([j|). 
It is not "universal" in the sense 
of describing any single halo drawn 
from a random cosmology (nor was 
it ever described as such by NFW, 
who point out its applicability to re- 
laxed systems). 

For all its newfound faults, all ex- 
periments done to date confirm that 
the form of equation(ij) does a re- 
markable job of describing a signifi- 
cant fraction of the mass density pro- 
files in the majority of collapsed ha- 
los. The NFW profile remains the 
de facto standard analytic form for 
halo density profiles, and it will re- 
main so until future studies provide 
well calibrated modifications. With 
apologies to Ivan King, a return to 
truncated isothermal spheres (King 
1962) would be a step in the wrong 
direction. 

One of the practical ramifications 
of Jing's work is that observers can- 
not rely on observation of the mass 
profile in a single cluster (e.g., Tyson, 
Kochanski & Dell'Antonio 1998) to 
place strong constraints on cosmol- 
ogy. Samples of ten or so clus- 
ters would provide interesting con- 
straints, particularly if the defining 
sample criteria were tailored to se- 
lect against objects undergoing sig- 
nificant mergers. X-ray imaging and 
spectroscopy could provide the clues 
necessary for such sample selection, 
but this requires understanding of 
the dynamical behavior of the intr- 
aclustcr medium in clusters. 



3 "Dissipationless" Gas Dy- 
namics : X— ray Clusters 

The hot intracluster plasma that 
permeates clusters of galaxies rep- 
resents the next most important 
matter constituent of clusters after 
dark matter. The fact that this in- 
tracluster medium (ICM) outweighs 
the matter associated with cluster 
galaxies is one of the legacies of 
the ROSAT satellite mission. Based 
on the ROSAT all-sky survey imag- 
ing of Briel, Henry & Bohringer 
(1992) and galaxy photometric data 
of Godwin, Metcalfe & Peach (1983), 
White etal. (1993) find a gas-to- 



galaxy mass ratio M gas /M ga i = 
(5.5 ± 1.5)?t- 3 '' 2 for the Coma clus- 
ter within an Abell radius ta = 
1.5 /i — 1 Mpc. This translates to a 
factor 10 for h = 0.65. 

The dominance of the ICM in the 
baryonic budget of clusters has a 
number of implications. From a mod- 
eling perspective, it supplies ammu- 
nition to justify the assumption that 
interaction between galaxies and the 
ICM is a "higher— order" effect, at 
least in very rich clusters such as 
Coma. This leads to the basic treat- 
ment of the dynamical and thermo- 
dynamic behavior of the ICM as be- 
ing driven by the gravitational evo- 
lution of the dominant, dark mat- 
ter. A complex network of shocks of 
varying strengths develop naturally 
in the merger /accretion process and 
the combined action of these shocks 
heats the gas by thermalizing the in- 
fall energy of the gravitational clus- 
tering process. 

In the nearly decade since this pro- 
cess was first examined using self- 
consistent collisional gas and col- 
lisionless dark simulations matter 
(Evrard 1990a, b), much has changed 
but much has stayed the same. One 
of the things that has stayed the 
same is efficiency of thermalization 
of the gas, traditionally expressed 
as /3 = o^ D /(kT/ [imp), the ratio 
of specific energies in dark matter 
and gas. In a comparison of twelve 
gas dynamics codes to the forma- 
tion of a single X-ray cluster (drawn 
from a standard cold dark matter 
model with 10% baryons), Frenk 
etal. (1998) find (fi) = 1.17 ± 0.06 
within r2oo to be the mean and 
standard deviation among the twelve 
codes, consistent with the value /3 = 
1.2 found by Evrard (1990a). A note- 
worthy point is that, despite their 
very different numerical treatments 
of shocks, codes of different type — 
Eulerian, Lagrangian hybrid — pro- 
duced very similar values of (3. 

A look at the P3MSPH code 
(Evrard 1988) solution to the com- 
parison cluster is provided in Fig- 
ures U and ^ which show four 
principle cluster observables at red- 
shifts z = 1.2, 0.6 and 0, respectively. 
The four panels show greyscale rep- 
resentations of the logarithms of 
the anticipated thermal Sunyaev— 
Zel'dovich signal (over 2 decades), 
intrinsic X-ray surface brightness 
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Figure 4. Views of a cluster at z = 1.24. 



(4 decades), gas temperature (fac- 
tor 4, from 0.25 - 1 X 10 8 K) and 
dark matter (DM) surface density 
(2.5 decades). The gas maps are 
smoothed on the local SPH ker- 
nel scale while the DM is left un- 
smoothcd to show the graininess in 
the particle solution. A physical re- 
gion of size 3.9 h -1 Mpc is shown at 
all redshifts. An mpeg animation of 
this four— panel figure is available a t 
http: / /astro. physics. lsa.umich.edu/ 
evrard/VCL/ccp_movie.htm. 

There are characteristic features 
in these images common to all twelve 
codes in the comparison study. At 
z = 1.2. most of the mass in the 
region has already collected into a 
long, knotty filament which will set 
the orientation for a chain of future 
merger events between the knots. 
One can imagine a network of such 
filaments covering the sky at these 
redshift — the "proto-cluster for- 
mation epoch". A good portion of 
the absolute scale the SZ decrement 
shown in the figures, 3.4 to 340 ^K, 
is visible with current ground based 
techniques (e.g., Carlstrom, Joy & 
Grego 1996) and the cosmic web that 
is the large-scale distribution of clus- 



ters (Bond & Myers 1998) may be re- 
vealed by future CMB experiments. 

An interesting set of indepen- 
dently observable features occurs in 
a pre— merger phase displayed in Fig- 
ure |i| The merger event involves a 
triplet, consisting of a pair of satel- 
lites (projected out along the line- 
of-sight in this image) each with 
mass in excess of 20% of the dom- 
inant precursor's mass. At the time 
shown, the thermalized gas envelopes 
of the triplet are being compressed 
and mildly shocked. The X-ray im- 
age displays two prominant peaks 
(the projected pair of satellites is 
to the upper right) and a strong 
emission weighted temperature gra- 
dient exists along the collision axis 
near its vertex. This pinching effect 
comes about because the gas at the 
cores of the precursors has not yet 
been moved off their incoming adia- 
bats, while gas trapped near the sys- 
tems center of mass is squeezed hard 
and subsequently strongly shocked. 
The hot pinched spot with adjacent 
sidelobes is a distinct signature of a 
merger occurring along the line of 
sight. 

The electron pressure visible in the 
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Figure 6. Views of a cluster at z = 0. 
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SZ image has features common to 
both the gas density (X-ray) and 
temperature. Two peaks are sepa- 
rated by a pressure ridge of strongly 
shocked material perpendicular to 
the collision axis. Because of the con- 
fining ram pressure along the col- 
lision axis, the flow of gas at this 
time is highly anisotropic. Inward 
flow along the axis feeds outward 
flow off-axis, and a small fraction of 
gas is shed into the voids which sur- 
round the large-scale filament. The 
dark matter density distribution is 
highly flattened, with two peaks evi- 
dent as in the X-ray image. 

Observations with the same dy- 
namic range and resolution of these 
maps do not yet exist. In this sense, 
theory is ahead of the observation, 
and one can interpret these maps as 
predictions for upcoming, high reso- 
lution ground-based and satellite ob- 
servations. The spectroscopic imag- 
ing capabilities offered by AXAF will 
soon begin probing the line-of-sight 
velocity structure, and observable ve- 
locity splittings are anticipated (Nor- 
man & Bryan 1998). 

Despite this seemingly active for- 
mation history, relatively little of 
the gas that was initially associ- 
ated with the dark matter in the 
virial region is lost from the sys- 
tem at the present epoch. Figure ^ 
shows the normalized gas fraction 
T(r) = M g as(r)/n b Mtot(r) evalu- 
ated within r200 for the twelve clus- 
ters simulated in the comparison 
study. The solid and dashed lines 
show the mean and la range (T) = 
0.92 ± 0.06 determined from the 
experiments. The Lagrangian codes 
(open symbols) seem to lose some- 
what more of the gas than the other 
code types, but the source of this dif- 
ference is not yet understood. 

Determination of the normaliza- 
tion of this effect is important for 
constraints on the total mass density 
fl m from the mean cluster baryon 
fraction (White et al. 1993; Evrard 

1997) . Another important factor in 
that argument is the normalization 
of the virial mass— temperature rela- 
tion. Experiments show small scat- 
ter about the virial relation (Evrard, 
Metzler & Navarro 1996; Bryan & 
Norman 1998; Eke, Navarro & Frenk 

1998) , but the zero point remains un- 
certain to perhaps 15%. A look at the 
agreement between virial mass and 
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Figure 7. Normalized gas fractions 
within r2oo f or the twelve codes 
in the cluster comparison study. 
Open symbols are Lagrangian (SPH) 
codes, filled circles Eulerian and 
stars hybrid. 
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Figure 8. Temperatures and masses 
from the cluster comparison study. 
Note the 20 — 30% range of the axes. 



(mass-weighted) temperature for the 
cluster comparison study is given in 
Figure H. Hints of the overall relation 
M oc T 5 / 3 in this diagram suggest 
that part of the disagreement among 
codes is due to a lack of strict syn- 
chronization in the cluster's dynam- 
ical state. Small differences in the 
treatment of tidal fields in the linear 
regime, for example, can be amplified 
into large orbital phase differences 
for infalling small satellites such as 
that visible to the north in Figure ^. 
Images based on the twelve codes' so- 
lutions exhibit noticeable differences 
in the positions of satellite objects. 

Table ^ summarizes the extent of 
the agreement among codes for inter- 
nal structural properties of the clus- 
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Table 2. Level of agreement between 
twelve different cosmological gas dy- 
namics codes in the final, virial 
structure of a single X-ray cluster 
(Frenk etal. 1998). 



Quantity 


rms deviation 


Pdm 


0.031 


Pgas 


0.086 


&DM 


0.046 


Tgas 


0.053 


Pgas 


0.081 


Lx 


0.23 



ter at the final time. The rms devi- 
ation of spherically averaged, binned 
data interior to the virial radius of 
the cluster (using 0.2 dex binwidth) 
is shown. Agreement in all the listed 
quantities is better than 10%, except 
for the bolometric X— ray luminosity, 
which has an rms of 23%. The lumi- 
nosity h x = J dVp 2 gaa (r)T^ 2 (r) is 
particularly sensitive to the core gas 
structure, and the larger deviation in 
Lx reflects the fact that the codes 
agree less well on the core structure 
than they do on the exterior enve- 
lope. 

Additional physics such as radia- 
tive cooling and magnetic fields will 
also affect the core structure, per- 
haps quite strongly. These effects are 
unlikely to affect the structure of the 
bulk of the gas in very rich clusters, 
but may contribute perhaps tens of 
percent effects. Many effects are only 
now being considered with careful 
experiments. Dolag, Bartelmann & 
Lcsch (1998) displayed a poster at 
the meeting showing that the effects 
of magnetic fields on cluster evolu- 
tion. Teyssier, Chieze & Alimi (1997) 
include non-equilibrium thermody- 
namics and examine differences in 
electron and ion temperatures. Met- 
zler & Evrard (1997) examine the ef- 
fects of galactic winds. None of this 
additional physics appears to change 
the overall picture of the ICM pre- 
sented above, though some details 
will certainly be affected. Dissenters 
point out that a strongly multi- 
phase intracluster medium is not for- 
mally ruled out (Gunn &; Thomas 
1996), but upcoming SZ maps com- 
bined with AXAF and XMM spec- 
troscopic imaging will provide defini- 



tive constraints (Nagai, Sulkanen & 
Evrard 1998). 

4 Galaxy Formation 

Galaxy formation is the perennial 
frontier in the simulation business. 
As outlined in Figure the radiative 
cooling instability within collapsed, 
dark matter potential wells is the es- 
sential ingredient which allows the 
baryons to sink to the center, be- 
come self— gravitating and, thereby, 
light up (White k, Rees 1978). Re- 
solving this instability is numerically 
challenging (see the heroic efforts of 
Abel, Bryan & Norman in this vol- 
ume). One has the added complica- 
tion that, even with a "perfect" nu- 
merical code, the physics governing 
star forming regions is not yet under- 
stood. The problem is formally ill- 
posed and solutions will be guided as 
much by instinct and intuition as by 
first principles. 

Radiative cooling of the gas has 
been enabled in a number of calcu- 
lations (Katz & Gunn 1991; Katz, 
Hcrnquist & Weinberg 1992; Evrard, 
Summers & Davis 1994; Navarro 
k. White 1994; Steinmetz & Miiller 
1994; Frenk etal. 1996; Tissera, 
Lambas & Abadi 1997; Navarro & 
Steinmetz 1997) and a population of 
cold, condensed baryonic cores de- 
velops within their enveloping halos, 
largely along the lines scripted by 
the sages twenty years ago. A num- 
ber of potentially interesting results 
have emerged, but their robustness 
remains in question. 

The issue of "biasing" — how 
the phase space structure of galax- 
ies differs from that of the dark 
matter — is of crucial importance 
to interpreting data from upcom- 
ing deep optical surveys such as the 
Two Degree Field (2dF, see Mad- 
dox in this volume) and the Sloan 
Digital Sky Survey (SDSS). Bias- 
ing is a complicated issue because 
both gravitational (dynamical fric- 
tion, merging, galaxy harassment) 
and non-gravitational effects (de- 
pendence of gas cooling and star for- 
mation rates on large— scale environ- 
ment) may play equally important 
roles. Direct simulation of cluster en- 
vironments {e.g., Frenk etal. 1996) 
in conjunction with semi-analytic 
approaches which postulate recipes 
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a galaxy in a diffuse dark mailer halo 



Figure 9. Flow diagram of the galaxy formation process. 



for star formation and feedback in 
halos (e.g., Kauffmann, Nusser & 
Steinmetez 1997) are beginning to 
yield clues. 

An example of a result which 
may turn out to be robust is shown 
in Figure |lo[ The figure shows the 
pairwise peculiar velocity of galax- 
ies and dark matter as a func- 
tion of pair separation derived from 
a 16 million particle simulation of 
an LCDM universe (Pearce et al. , 
in preparation). Galaxies — identi- 
fied with cold, high density baryonic 
knots in the calculation — show a 
lowered overall velocity distribution 
compared to the dark matter and 
this produces a much better agree- 
ment with observational values de- 
rived from the LCRS redshift sur- 
vey (Lin etal. 1996). Whether or not 
this agreement is fortuitous hinges on 
the sensitivity of this result to details 
such as numerical parameter choices, 
galaxy sample definition (both real 
and simulated) and the weighting of 
the velocity statistic (Davis, Miller & 
White 1997). Still, the agreement be- 
tween theory and observation in Fig- 
ure is encouraging and it may in- 
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Figure 10. Distribution of peculiar 
velocities for galaxies derived from a 
16 million particle SPH calculation 
(Pearce etal. , in preparation). 

dicate that the current simulations 
are accurately capturing galaxy dy- 
namics within the large— scale web of 
dominant dark matter. 
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5 Summary 

The present state of numerical 
modeling of large— scale structure in 
the universe is 

• well understood and ahead of ob- 
servations of dark matter structure, 

• becoming understood and on par 
with observations of X— ray clusters 
and the Ly-a forest, 

• poorly understood and far be- 
hind observations of the galaxy dis- 
tribution. 

What does the future hold? We 
can expect to see parallel computers 
used both to expand dynamic range 
(by allowing bigger calculations) and 
also to explore model parameter 
space (by allowing large numbers of 
smaller calculations). The number 
and range of physical processes being 
included in the calculations will cer- 
tainly increase. The downside of this 
added complexity is that it will nat- 
urally introduce additional sources 
for numerical error. Besides highly 
simplified systems, the main check 
on the numerical solutions may re- 
main internal — a "galaxy compar- 
ison project" analogous to that per- 
formed for clusters is likely to occur 
in the near future. The upside of the 
added complexity is the improved 
contact it will bring with observa- 
tional data. Spectroscopic imaging, 
broad band colors, kinematic studies, 
etc. — the sky is the target (but not 
necessarily the limit) in the virtual 
world. 

It is easy to foresee the prob- 
lem of galaxy formation being scru- 
tinized more thoroughly both from 
the "bottom— up" (via modeling of 
absorption line systems and high-z 
progenitors of galaxies) and from the 
"top-down" (via modeling of galax- 
ies in different large-scale environ- 
ments, clusters versus the field). By 
attacking the problem of galaxy for- 
mation from below and above in this 
way, we will squeeze out the answer 
to the lingering question "how do 
galaxies trace the mass" in the uni- 
verse. 
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